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We present a detailed quantum mechanical non empirical DFT investigation of the energy- 
optimized geometries, phase stabilities and electronic properties of bulk PtaN, PtN and PtN2 in 
a set of twenty different crystal structures. Structural preferences for these three stoichiometries 
were analyzed and equilibrium structural parameters were determined. We carefully investigated 
the band-structure and density of states of the relatively most stable phases. Further, GWo cal- 
culations within the random-phase approximation (RPA) to the dielectric tensor were carried out 
to derive their frequency-dependent optical constants of the most likely candidates for the true 
crystal structure. Obtained results were comprehensively compared to previous calculations and to 
experimental data. 
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I. INTRODUCTION 

Platinum is known to form simple binary compounds 
with other elements (e.g. PtF 4 , Ptl 2 , PtO and PtS) 1 . 
However, platinum had not been known to form crys- 
talline solid nitride, but other forms of platinum nitrides 
(e.g. PtN 1 - 2 , PtN 2 \ (PtNh 1 ' 2 and Pt 2 N 2 ) had been ob- 
served. 

In January 2004, Soto 3 reported the preparation of 
platinum thin films containing up to ~ 14 at. % nitro- 
gen. The study concluded that platinum can form an 
incipient nitride phase with composition near to Pt@N. 
Few months later, in May 2004, Gregoryanz and co- 
workers 1 published the discovery and characterization of 



solid crystalline platinum mono-nitride for the first time. 
The synthesis was carried out above 45 GPa and 2000 K 
but with complete recovery of the product at room pres- 
sure and temperature. The produced samples have a very 
high bulk modulus leading to important implications in 
high-pressure physics and technology. The 1 : 1 stoichiom- 
etry was assigned to the new nitride, and according to 
their XRD measurements, Gregoryanz et al. proposed 
three structures: Bl, B3 and B17 (for description of the 
structures see sub-section II A below) , all based on the Pt 
fee sub-lattice. Due to some considerations, Bl and B17 
were ruled out and B3 was assigned to the new product. 

In addition to the well-crystallized and highly or- 
dered regions, a common feature in the synthesized 
platinum nitrides is the presence of sub- or/and super- 
stoichiomctric phases containing N or Pt vacancies and 
residual non-stoichiomctric material distributed through- 
out the samples 1,4 . 

The work of Gregoryanz et al. 1 has stimulated many 
further theoretical studies 5-7 as expected by Gregoryanz 
and co-workers 1 themselves. However, theoretical work 
showed that PtN(B3) is elastically and thermodynami- 
cally unstable (see sub-section III A below). Accordingly, 
claiming that large errors are generally inevitable in the 
used experimental characterization methods 6 ' 8 , and due 
to other paradoxical facts 9 in the original experiment by 
Gregoryanz et al. 1 , theoreticians questioned the chem- 
ical stoichiometry and the crystal structure of this new 
material and started to investigate other possibilities 8-10 . 
Moreover, the experimentally reported 1 high bulk mod- 
ulus of the platinum nitride has not been reproduced by 
any reliable calculations and its mechanism is still an un- 
clear open problem 5 ' 9,10 . 

These investigations led to a kind of consensus that the 
compound does not crystallize in the proposed PtN(B3) 
phase 4 , but the true stoichiometry and the true crystal 
structure have become now a matter of debate 6 ' 9 . 

In an apparent attempt to respond to this debate, 
Crowhurst et al. 4 managed, in 2006, to to reproduce and 
characterize platinum nitride. Combining theory with 
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their own observed Raman spectrum, they came up with 
a conclusion to propose PtN 2 (C2) and rejected PtN(B3), 
proposed by the first platinum nitride synthesizers 1 , and 
PtN2(Cl), proposed in some theoretical works. Like the 
first proposed structures 1 , CI 8 and C2 4 structures have 
the fee sub-lattice of the metal. 

Despite the considerable number of the subsequent the- 
oretical studies, the discrepancy between theory and ex- 
periment in the structural and the physical properties 
of this nitride is not yet satisfactorily understood. Nev- 
ertheless, many transition metals can form more than 
one nitride 11 . Thus, it is of interest to know if platinum 
can form nitrides with different stoichiometrics and/or 
crystal structures other than those proposed by the first 
platinum nitride synthesizers and other researchers. 

In the present work, we present a comprehensive first- 
principles calculations of the equation of state, possible 
pressure-induced phase transitions, electronic and opti- 
cal properties of crystalline PtsN, PtN and PtN 2 in a 
total of twenty different -previously proposed and new 
hypothetical- structural modifications. The work partly 
aims to solve some of the reported discrepancies. In ad- 
dition, to the best of our knowledge, there is no available 
experimental or theoretical optical data for the platinum 
nitride, and the present study may be the first one to 
calculate the optical spectra of platinum nitrides. 



II. CALCULATION METHODS 
A. Stoichiometries and Crystal Structures 

The structure and stability of solids arc influenced by 
their chemical stoichiometry 12 and the electronic struc- 
ture of the outer shells of atoms is a controlling factor in 
proposing any crystal structure 13 . Like other theoretical 
works (cf. Table I) we postulated various structure types 
that are not based on the observed fee Pt sub-lattice. 
Our assignment of the following different chemical sto- 
ichiometries and crystal structures is based on the fact 
that many transition-metal nitrides (TMNs) are known 
to form more than one nitride 11 . Ni, which shares the 
same group with Pt, and Au which shares the same pe- 
riod with Pt in the periodic table, are known to form 
N13N and AU3N nitrides. Thus, it is of interest to con- 
sider more Pt atoms in the unit cell and less symmet- 
ric structures, and to know whether platinum can form 
PtsN with the reported structures of these and other 3:1 
TMNs. In this work, the hypothetical PtaN is studied in 
the following seven structures: DO3 (space group Fm3m), 
A15 (space group Pm3n), DO9 (space group Pm3m), Ll 2 
(space group Pm3m), D0 2 (space group Im3), e-FeaN 
(space group P6322) and RI1F3 (space group R3c). 

The following nine structures were assigned to PtN: Bl 
(space group Fm3m), B2 (space group Pm3m), B3 (space 
group F43m), B81 (space group P63/mmc), Bk (space 
group P63/mmc), B^ (space group P6m2), B4 (space 
group P63ULC), B17 (space group P4 2 /mmc) and B24 



(space group Fmmm). While CI (space group Fm3m), 
C2 (space group Pa3), C18 (space group Pnnm) and 
CoSb 2 (space group P2i/c) are the four structures which 
were proposed for PtN 2 . 



B. Electronic Relaxation Details 

As implemented in the all-electron Vienna ab initio 
Simulation Package (VASP) 14-19 , our electronic struc- 
ture calculations were based on spin density functional 
theory (SDFT) 20 ' 21 . To solve the self-consistent Kohn- 
Sham (KS) equations 22 

where i, k and a are the band, k-point and spin indices, 
respectively, the pseudo part of the KS orbitals ^ ,k (r) 
are expanded on plane- waves (PWs) basis. Only those 
PWs with kinetic energy ^-|k + G| < E cut = 600 eV 
were included. This always corresponds to changes in the 
total electronic energy and in the- so-called Fermi energy 
Ep that are less than 3 me V/ atom and 1 m.eV, respec- 
tively. 

For ionic relaxation, the Brillouin zones (BZs) were 
sampled using 13x13x13 T-centered Monkhorst-Pack 
meshes 23 ; while for the static calculations, 17 x 17 x 17 
meshes were used. Any increase in the density of this 
mesh corresponds to change in the total energy and in 
Fermi energy Ep that are less than 2 me^/atom and 
0.02 eV, respectively. 

Partial occupancies were set using the tetrahedron 
method with Blochl corrections 24-26 for the static total 
energy and the electronic density of states (DOS) calcula- 
tions; while in the ionic relaxation, the smearing method 
of Methfessel-Paxton (MP) 27 was used. In the latter case, 
the smearing width was set such that the fictitious en- 
tropy is always less than 1 me V/ atom. 

The Perdew, Burke and Ernzerhof (PBE) 28 ~ 30 
parametrization of the generalized gradient approxi- 
mation (GGA) 31 ~ 33 was employed for the exchange- 
correlation potentials V^-' c [n(r)}. VASP treats the core- 
valence interactions, V ex t(r) using the projector aug- 
mented wave (PAW) method 19 ' 34 . The PAW potential 
explicitly treats the 2s 2 2p 3 electrons of N and the ScPGs 1 
electrons of Pt as valence electrons. 

In the standard mode of VASP, while fully relativistic 
calculations are performed for the core-electrons, only 
scalar kinematic effects are incorporated to treat the 
valance electronic structure 18 . It was found that this 
scheme is sufficient and the spin-orbit interactions have 
little effect on the macro-physical properties of platinum 
nitride 5 . Thus, we made no effort to consider spin-orbit 
coupling of the valence electrons. 
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The relaxation of the electronic degrees of freedom was 
performed using the blocked Davidson iteration scheme 35 
as implemented in VASP. The electronic self-consistent 
(SC) convergence was considered to be achieved when the 
difference in the eigenvalues and in the total energy be- 
tween two successive steps are both less than 1 x 10~ 4 eV. 



C. Geometry Optimization and EOS 

At a set of isotropically varying volumes of the unit 
cells, ions with internal free parameters were relaxed un- 
til all Hellmann-Feynman force components 36 on each ion 
were less than 1 x 10~ 2 eV/K. Static total energy calcu- 
lations (as described above) at each volume followed and 
the cohesive energy per atom was calculated from 37,38 
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Here, Z is the number of Pt m N„ per unit cell, E^ om and 
^atom are the energies of the spin-polarized non-spherical 
isolated Pt and N atoms, E^^ n are the bulk cohesive 
energies calculated by VASP with respect to spherical 
non spin-polarized reference atoms, and m, n = 1, 2 or 3 
are the stoichiometric weights. 

The calculated E co h per atom as a function of volume 
V per atom were fitted to a Birch-Murnaghan 3rd-order 
equation of state (EOS) 39 



E(V) = E Q + 
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The four equilibrium fitting parameters: the equilibrium 
volume Vb , the equilibrium cohesive energy Eq , the equi- 
librium bulk modulus 
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and its pressure derivative 
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were determined by a least-squares method. 



D. Formation Energy 



Ef can be obtained from 



solid 



= i?coh(Pt m N^ lid ) 

m£ coll (Pt solid ) + 



f£ coh (Nf s ) 



(7) 



where m,n = 1,2,3 are the stoichiometric weights and 
-Eco^PtmN^ 01 " 1 ) is the cohesive energy per atom as ob- 
tained from Eq. 2. The ground-state cohesive en- 
ergy and other equilibrium properties of the elemen- 
tal platinum E^ (Pt solld ) in its fee Al structure (space 
group Fm3m No. 225) 40 ~ 42 are given in Table I. We 
found the equilibrium cohesive energy of the molecular 
nitrogen (i? co h(N| as )) and its N-N bond length to be 
—5.196 eV/atom and 1.113 A. For details on how these 
properties were calculated, readers are referred to Ref. 
43. 



E. GWA Calculations and Optical Properties 

Accurate quantitative description of optical proper- 
ties of materials requires treatments beyond the level of 
DFT 44 . One choice is to follow the method which is pro- 
vided by the many-body perturbation theory (MBPT). 
In this approach one needs to solve a system of quasi- 
particlc (QP) equations 45-47 

In practice, one takes the wave functions ip?^(v) from 
the DFT calculations. However, this technique is com- 
putationally expensive, and we had to use less dense 
meshes of k-points: (10 x 10 x 10) in the case of B17 and 
(12 x 12 x 12) in the case of B3. The quantity £(r, r'; ef£) 
in Eqs. 8 above is known as self-energy. It contains all 
the static and dynamic exchange and correlation effects, 
including those neglected at the DFT-GGA level. When 
£ is written in terms of the Green's function G and the 
frequency-dependent screened Coulomb interaction W as 



J de'G(r,r';e,e')W(ry; 



(9) 



the approximation is referred to as GW approximation. 
The dynamically screened interaction W is related to the 
bare interaction v via 



In addition to E co \ n another important measure of rel- 
ative stability is the so-called formation energy Ef. As- 
suming that the solid Pt m N„ results from the interac- 
tion between the solid Pt(/cc) metal and the gaseous N2 
through the chemical reaction 

>r ^ pt m Nf id , (6) 



mPt s 



Dlici 



W(r,r';e)=j / dne^r, n; e)v(n, r'), (10) 



where the dielectric matrix e is calculated within the ran- 
dom phase approximation (RPA). The QP eigenvalues 



Re 



OP 
k 



Hks — Vxc + ^GW 
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are updated in the calculations of G, while W is kept 
at the DFT-RPA level. This is called the GW self- 
consistent routine on G. After the execution of the fourth 
iteration, e is recalculated within the RPA using the up- 
dated QP eigenvalues 46-48 . It is straightforward then to 
calculate all the frequency-dependent optical spectra (e. 
g. refractive index n (w), extinction coefficient k (uj), ab- 
sorption coefficient a (uj), reflectivity R (uj) and transmi- 
tivity T (uj) = 1 — R(uj)) from the real £ Te (ui) and the 
imaginary £j m (u;) parts of £rpa(w) 49 ~ 51 : 



n (uj) 
k (uj) 
a (uj) 
R(uj) 



^| ( [4 (w) + e? m (uj)] » + e rc (w)) ' (12) 
i([e? c H + iH]^e rc H) 5 (13) 



[e re (Uj) + j£im (uj)} 2 - 1 



[e re (uj) + je im (w)p 



(14) 
(15) 



It may be worth to emphasize here that for more ac- 
curate optical properties (e.g. more accurate amplitudes 
and positions of the characteristic peaks), electron- hole 
excitations should be calculated by solving the so-called 
Bethe-Salpeter equation, the equation of motion of the 
two-body Green function G%- The latter can be evalu- 
ated on the basis of our obtained GW one-particle Green 
function G and QP energies 52 . 



III. RESULTS AND DISCUSSION 

Cohesive energy E CQ h versus atomic volume Vq equa- 
tion of state (EOS) for the different phases of Pt 3 N, PtN 2 
and PtN are displayed graphically in Fig. 1, Fig. 2 and 
Fig. 3, respectively. The corresponding obtained equi- 
librium structural parameters and energetic and clastic 
properties are presented in Table I. In this table, as well 
as in Fig. 4, structures are first grouped according to the 
nitrogen content, starting with the stoichiomctry with 
the lowest nitrogen content PtaN, followed by the 1:1 
series and ending with the nitrogen-richest PtN 2 group. 
Within each series, structures are ordered according to 
their structural symmetry, starting from the highest sym- 
metry (i.e. the highest space group number) to the least 
symmetry. Whenever possible, our results are compared 
with experiment and with previous calculations. In the 
latter case, the calculations methods and the XC func- 
tionals are indicated in the Table footnotes. 

To study the effect of nitridation on the elemental 
Pt(Al) and to easly compare the properties of these 
phases relative to each other, the calculated equilibrium 
properties are displayed relative to the corresponding 
ones of Pt(Al) in Fig. 4. 
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FIG. 1. (Color online.) Cohesive energy _E co h(el//atom) ver- 
sus atomic volume V (A 3 /atom) for PtaN in seven different 
structural phases. 
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FIG. 2. (Color online.) Cohesive energy _E co h(eV/atom) ver- 
sus atomic volume V (A 3 /atom) for PtN in nine different 
structural phases. 




FIG. 3. (Color online.) Cohesive energy _E co f l (e 1 F/atom) ver- 
sus atomic volume V(A 3 /atom) for PtN2 in four different 
structural phases. 
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TABLE I. The calculated (Pres.) zero-pressure properties of the platinum nitrides and the results of previous calculations 
(Comp.). The available experimental data are given in the last row. 







a(A) 


6(A) 


c(A) 


a(°) or Vo(A 3 /atom) B coh (eV/atom) B (GPa) 


Bo 


Ef (eVyatom) 


Pt 




Pres. 


3.978 






15.74 


-5.451 


242.999 


5.486 






Expt. 


(3.9233 ± 0.0007)" 






15.097° 


-5.84° 


278.3°, 280 d 


5.18° 




Al 


3.924"° 










249°° 


5.23" 






Comp. 


3.90 s , 3.890', 3.981 s8 , 








-7.04", 


-3.74', 305 s , 320 s , 238* 


5.16 k , 5.30 1 , 5.25 m 






3.967', 3.966" 








-5.53' 


249", 242 dd , 244.18 s8 


5.23", 5.83 dd , 5.7 SS 




Pt^N 


D0 3 


Pres. 


6.106 






14.23 


-4.140 


218.097 


5.282 


1.247 


A15 


Pres. 


4.924 






14.92 


-3.759 


194.136 


5.266 


1.628 


D0 9 


Pres. 


4.114 






17.41 


-4.558 


167.839 


5.241 


0.829 


Ll 2 


Pres. 


3.863 






14.41 


-4.021 


205.279 


5.472 


1.366 


D0 2 


Pres. 


7.875 






15.26 


-4.644 


117.174 


12.098 


0.743 


c-Fe 3 N 


Pres. 


5.680 




5.293 


18.49 


-4.713 


217.035 


6.779 


0.674 


RhF 3 


Pres. 


5.463 






a = 58.640 13.97 


-4.688 


224.419 


5.412 


0.699 


PtN 




Pres. 


4.495 






11.35 


-3.945 


230.869 


5.059 


1.378 


Bl 




4.45 x , 4.50\ 4.41 z 










232 x , 230", 288 z 






Comp. 


4.471" 
4.491 ss 






10.66" 




251", 242 dd , 294 ff 
229.76 ss 


4.00", 4.78 dd 
4.9 SS 


1.365" 

[0.375 + E f (BlT)] hh 


B2 


Pres. 


2.819 






11.20 


-3.522 


238.187 


5.070 


1.801 


Comp. 


2.818 ss 










234.88 s8 


5.1 ss 






Pres. 


4.782 






13.67 


-4.203 


193.466 


5.031 


1.120 






4.7217°, 4.8250°, 4.779 ss , 










243.3°, 196.3°, 190.61 s8 


5.1 ss , 




B3 


Comp. 


4.6833°, 4.7889°, 4.8114 r , 










271.9°, 192.7°, 184', 




0.95" 




4.692 s , 4.780*, 4.760", 










244 s , 194', 213", 


4.00°° 


1.1" 






4.80 x ' y , 4.70 z , 4.699 hh 










192 x , 190 y , 232 z , 217 dd 


3.62 dd 


[0.21 + E f (Bn)] hh 


B81 


Pres. 


3.482 




4.843 


12.71 


-3.713 


210.165 


4.945 


1.610 


B k 


Pres. 


3.378 




8.986 


22.20 


-4.061 


108.968 


4.553 


1.262 


B h 


Pres. 


3.039 




2.966 


11.86 


-3.716 


222.279 


5.014 


1.607 


B4 


Pres. 


3.382 




5.539 


13.72 


-4.171 


190.130 


5.033 


1.152 


Comp. 


3.386 ss 




5.529 ss 






191.00 s8 


4.7 SS 




B17 


Pres. 


3.069 




5.403 


12.72 


-4.652 


235.041 


5.018 


0.671 


Comp. 


3.323 hh 




4.579 hh 












B24 


Pres. 


4.216 


4.472 


4.948 


11.66 


-3.928 


226.608 


5.153 


1.395 


Comp.* 3.972 hh 


3.977°° 


6.022 hh 






270 hh 




[0.085 + E J (B17)\ hh 


PtN 2 




Pres. 


4.963 






10.19 


-3.918 


263.295 


4.717 


1.363 


CI 


Comp. 


4.9428°, 5.0403°, 










322.1°, 267.2°, 269", 


4.00°° 


1.167°" 




4.866 s , 4.958*, 4.939 dd 










316", 264', 260 dd 


4.73 dd 


1.317 d ° 




Pres. 


4.912 






9.882 


-4.689 


226.779 


6.893 


0.592 


C2 


Comp. 


4.87" 






9.12 H 








0.267", 0.24 bb 




4.848", 4.874" 






9.65" 




305", 285 dd , 300" 


4.00", 5.50 dd 


0.64", 0.212" 


C18 


Pres. 


3.036 


3.984 


4.862 


9.800 


-4.755 


244.320 


7.938 


0.526 


Comp. 


3.778" 


4.880" 


3.208" 


9.827" 




286" 




0.249" 


CoSb 2 


Pres. 


5.460 


5.163 


9.374 


13 = 151.225 10.60 


-4.508 


118.594 


6.619 


0.773 


Comp. 


4.950" 


4.880" 


4.950" 


99.50" 9.827" 




289" 




0.248" 


Experiment 






i 1.NIM2 X 5)"- v . 1.80 11(2)" 










(372 ±5)", (354 ±5) v 


4.0°, 5.26" 





a Ref. 53: This is an average of 23 experimental values, at room temperature. 
b Ref. 53: at room temperature. 

c Ref. 54: Cohesive energies are given at A" and 1 atm = 0.00010 GPa; while bulk moduli are given at room temperature. 
d Ref. (25) in 55: at room temperature. 
e See Refs. (8)-(ll) in 55. 

g Ref. 56: using the full-potential linearized augmented plane waves (LAPW) method within LDA. 
h Ref. 57: using the projector augmented wave (PAW) method within LDA. 
1 Ref. 57: using the projector augmented wave (PAW) method within GGA(PW91). 
J Ref. 57: using the projector augmented wave (PAW) method within GGA(PBE). 

k Ref. 55: using the so-called method of transition metal pscudopotential theory; a modified form of a method proposed by Wills and 

Harrison to represent the effective interatomic interaction. 
1 Ref. 55: using a semi-empirical estimate based on the calculation of the slope of the shock velocity vs. particle velocity curves 
obtained from the dynamic high-pressure experiments. The given values are estimated at ~ 298 K. 
m Ref. 55: using a semi-empirical method in which the experimental static P — V data are fitted to an EOS form. The given values are 
estimated at ~ 298 K. 

n Ref. 5: using the ultrasoft pseudopotcntial (USPP) method within LDA. Bo's are calculated from elastic constants. 

Ref. 5: using the ultrasoft pseudopotcntial (USPP) method within GGA. Bo's are calculated from elastic constants. 

p Ref. 5: using the projector augmented wave (PAW) method within LDA. B' is set to be 4. 

q Ref. 5: using the projector augmented wave (PAW) method within GGA. B' is set to be 4. 

r Ref. 5: using fully relativistic full-potential linearized augmented plane waves (LAPW) method within GGA. 

s Ref. 8: using the full-potential linearized augmented plane waves (LAPW) method within LDA. 

* Ref. 8: using the full-potential linearized augmented plane waves (LAPW) method within GGA(PBE). 

u Ref. 1: The experimental evolution of the volume with pressure was fitted with a Birch-Murnaghan EOS, but BL was set to be 4. 

v Ref. 1: The experimental evolution of the volume with pressure was fitted with a Birch-Murnaghan EOS, but B Q pt = 5.26 was fixed. 

w Ref. 1: From XRD measurements at 0.1 MPa. 

x Ref. 9: using the full-potential linearized augmented plane waves (FPLAPW) method within GGA(PBE). 

y Ref. 9: using pseudopotentials method within GGA(PBE). 

z Ref. 9: using pseudopotentials method within LDA. 
aa Ref. 4: using the PAW method within GGA(PW91), but the experimental value of B C oh(N| as ) in Eq. 7 was used. 
bb Ref. 4: using the full-potential linear-augmented plane-wave method. 
cc Ref. 6: using pseudopotentials method within GGA(PBE). 
dd Ref. 6: using pseudopotentials method within GGA(PBE). 
ee Ref. 58. 

ff Ref. 7: using the pseudopotcntial method within LDA. 
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FIG. 4. (Color online.) Calculated equilibrium properties of the twenty studied phases of platinum nitrides. All quantities are 
given relative to the corresponding ones of the fee crystalline elemental platinum given in the first row of Table I. 
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A. EOS and Relative Stabilities 

Fig. 1 reveals that Pt 3 N in its least symmetric phase, 
the trigonal (rhombohedric) structure of RI1F3, is the 
most favorable phase in this series. However, after 

- 15.9 A 3 /atom the EOS of Pt 3 N(RhF 3 ) is almost iden- 
tical with the EOS of Pt 3 N in the simple cubic struc- 
ture of the anti-Re0 3 (DOg). They share a minimum at 
~ (17.4 A /atom, —4.56 eV). Very close to this point, 

at - (17.23 A 3 /atom, -4.59 eV), the EOS of Pt 3 N(D0 2 ) 
has a kink due to a change in the positions of some Pt 
ions. 

The EOS of Pt 3 N(Fc 3 N) has two minima 
located at (14.11 A 3 /atom, -4.697 eV) and 

(18.26 A 3 /atom, —4.679 eV). Thus, the two minima are 
very close in energy but, due to the difFcrcncc in Vo, they 
correspond to bulk moduli of 222.7 GP and 169.0 GP, 
respectively The Pt ions are in the 6g Wyckoff positions: 
(x, 0, 0), (0, x, 0), (-x, -x, 0), (-x, 0, |), (0, -x, |) and 
(x, x, 5). Upon ion relaxation of Pt 3 N(Fc 3 N), atomic 
positions change from x ~ 3 to x = \ causing the sudden 

change in the potential surface (at ~ 16.83 A 3 /atom) 
as the bulk Pt 3 N(Fe 3 N) being decompressed (Fig. 1). 
It may be worth mentioning here that Ag 3 N(Fe 3 N) 61 
and Cu 3 N(Fe 3 N) 43 were found to behave in a similar 
manner. 

Hence, one of the two minima in the EOS of 
Pt 3 N(RhF 3 ) is shared with the minimum of the EOS of 
Pt 3 N(D0g) and the other is shared with one of the two 
minima of Pt 3 N(Fe 3 N). 

The crossings of the less stable D0 3 , LI2 and A15 EOS 
curves with the more stable DO9 EOS at the left side 
of their equilibria indicates that DOg would not survive 
under pressure and that possible pressure-induces phase 
transitions from the latter phase to the former ones may 
occur. 

Fig. 4 shows that the Pt 3 N most stable phases may 
energetically compete with the PtN and PtN 2 most 
stable ones. However, from the foregoing discussion, 
it seems that Pt 3 N would not have a simple potential 
surface. 

Using the full potential augmented plane wave plus lo- 
cal orbitals (APW+lo) method within GGA(PBE), the 
energy-volume EOS's for Bl, B2, B3 and B4 have been 
studied by the authors of Ref. 59. Some of their obtained 
equilibrium properties are included and referred to in Ta- 
ble I. Within the considered parameter sub-space, our 
obtained EOS's (Fig. 2), relative stabilities, and equi- 
librium structural parameters and mechanical properties 
(Table I) are in excellent agreement with their findings. 
However, relaxing the c/a parameter, they obtained an 
additional EOS which lies below all the other considered 
ones, but its equilibrium Bq is significantly smaller. 

From Fig. 2, it is evident that PtN(B17) is the ener- 
getically most stable phase in the PtN series. The differ- 



ence in the equilibrium i? co h between PtN (B 17) and the 
next (less) stable phase, PtN(B3), is about 0.5 eV (Table 
I). This difference was found by other researchers 9 to be 
0.9-1.05 eV. The crossings of the EOS curve of B17 with 
some of those of less stability at the left side of their equi- 
libria reveals possible pressure-induces phase transitions. 
To closely investigate these transitions, we plot the cor- 
responding relations between enthalpy H = E(V) + PV 
and the imposed external pressure P. Possible transi- 
tions and the pressures at which they occur are carefully 
depicted. A point where two H(P) curves (of two modifi- 
cations with the same chemical stoichiometry 7 ) meet rep- 
resents a phase transition from the phase with the higher 
H to the one with the lower H 37 . From the H(P) dia- 
grams (not shown here) we found that PtN(B17) would 
transform to PtN(Bl, B2, B h or B24) at - 93 GPa, 
- 143 GPa, - 193 GPa or - 123 GPa, respectively. 

It may be worth to mention here a few points about 
this B17 structure: (i) It was theoretically predicted to be 
the ground-state structure of CuN 43 , AgN 61 , AuN 62 and 
PdN 63 . (ii) The same foregoing phase PtN-PtN struc- 
tural pressure-induced transitions have been predicted 
for PdN, but at relatively smaller pressures in the range 
(25.8 - 62.1 GPa) 63 . (iii) B17 is the structure of PtS 64 
and PtO 9 . (iv) It was found by other authors to be a 
possible ground state for PtN 9 . (v) The B17 structure 
has an fee Pt sub-lattice (as the synthesized platinum ni- 
tride), but it is tetragonal and the sub-lattice are highly 
distorted {c/a ^ y/3 versus c/a = \/2 for ideal fee), and 
probably because of this distortion it was rejected by the 
platinum nitride synthesizers 1 , (vi) Fig. 5 that B17 is 
energetically favorable over Bl and B3 at all pressures. 

Nevertheless, PtN(B17) was found to be elastically 
unstable 10 . 

Assuming 1:1 stoichiometry, the first platinum ni- 
tride synthesizers assigned the B3 structure for their 
product 1 . However, it was shown in the same work that 
PtN(B3) should break down or transform at pressures 
above 12 GPa. In agreement with this experimental 
prediction, Fig. 5 shows that PtN(B3) would not sur- 
vive at pressures above 19 GPa where the B3— >B1 phase 
transition occurs. Other theoretical works also predicted 
that Bl becomes more favorable than B3 structure above 
13.3 GPa 9 , ~ 15 GPa 7 , 16.5 GPa 9 , and 17.6 GPa 9 . 

Therefore, we support Ref. 9 on the judgment that, 
unless the PtN(B3) was formed upon depressurization, 
its production at 45 — 50 GPa 1 is questioned. Further, 
first-principles calculations showed that PtN(B3) is 
elastically unstable 5,6,8 , and that it may distort sponta- 
neously to a tetragonal lattice to lower the energy 8 . 

In the PtN 2 series, we can see from Table I and from 
Fig. 4 that PtN 2 in the simple orthorhombic structure 
of FeS 2 marcasite (C18) is the most stable phase, while 
the face-centered cubic structure of CaF 2 fluorite (CI) 
is significantly the least favorable structure. Yet, Fig. 5 
reveals that the latter PtN 2 (Cl) is more favorable than 
the proposed PtN(Bl, B3 and B17) at pressures above 
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FIG. 5. (Color online.) Enthalpy H vs. pressure P equation 
of state (EOS) for the most favorable PtaN phase (RI1F3), the 
three proposed PtN modifications (Bl, B3 and B17), and the 
four considered PtN2 structures in the present work (CI, C2, 
C18 and C0SD2). The arrows indicate the pressures at which 
curves cross each other. 

4 GPa, 14 GPa and 52 GPa, respectively. Others 9 found 
PtN2(Cl) to be more favorable than PtN(B3) at pres- 
sures above 30 GPa. 

In contrast to our enthalpy-pressure EOS's in Fig. 
5, Chen, Tse and Jiang 60 obtained an H(P) curve 
for C18 which lies always above the curve for C2 and 
coincides with the one of CoSb2. They concluded that 
C2 is the most stable structure among these three 
modifications. While we sticked to the original C18 
relative dimensions, it seems that Chen, Tse and Jiang 
tried to optimize the lattice parameters ratios (see Table 
I). However, the c : a : b ratio they obtained is very 
close to our a : b : c ratio, and the difference in Vq is 
less than 0.03 A /atom) 65 . Another difference is the 
atomic electronic configuration of Pt 5g? 8 6s 2 they used. 
Nevertheless, they agreed with us that in the — 60 GPa 
pressure range, no transition between these three phases 
occurs. 

Comparing the relative stability of the three most sta- 
ble compositions, we find from Table I and from Fig. 
4 that PtN2(C18) is the most favorable, followed by 
Pt 3 N(RhF 3 ), and the least stable phase is PtN(B17). 
However, the differences in their equilibrium E co h lies 
within a narrow range of 0.036 eV. Relative to their par- 
ent metal, all phases have higher -E co h, i.e. they are less 
bound than Pt(Al). Hence, we found, as other theoret- 
ical works 8 , that platinum nitride can be stabilized in 
stoichiomctries and structures other than that proposed 
by the first synthesizers 1 . 

In Ref. 9, the energy-volume EOS for Bl, B3, B17, 
CI, and C2 have been studied using DFT-GGA. Within 
this parameter sub-space, our obtained EOS's (Figs. 2, 
3 and 4) are in excellent agreement with the findings of 
9. From the relative enthalpy-pressure diagrams 66 , Ref. 
9 arrived at an astonishing result: the experimentally 



proposed PtN(B3) is an entirely unstable structure at 
any pressure. 

To closely study the non-zero pressure stoichiometric 
and structural preferences, we displayed in Fig. 5 the 
enthalpy H vs. pressure P equation of states (EOS) for 
the most favorable Pt 3 N phase (RhF 3 ), the three pre- 
viously proposed PtN modifications (Bl, B3 and B17), 
and the four considered PtN2 structures in the present 
work (CI, C2, C18 and CoSb2). The arrows indicate 
the pressures at which curves cross each other. From 
these curves, it is clear that PtN 2 (C18), followed by 
PtN2(C2), are the most energetically favorable phases at 
all pressure. At pressures above 10 GPa, PtN 2 (CoSb2) 
has lower enthalpy than the rest of the modifications, in- 
cluding PtN(B17) and PtN 2 (Cl). At pressures higher 
than 8 GPa, PtN(B17) becomes more favorable than 
Pt 3 N(RhF 3 ), but the former never competes behind 
52 GPa when PtN2(Cl) becomes more favorable. How- 
ever, Pt 3 N(RhF 3 ) is more stable than PtN(B3) at all 
pressures. In summary, Fig. 5 reveals that even if a PtN 
phase has been observed (at pressures around 50 GPa), 
this phase must be unstable toward phase decomposition 
into solid constituents Pt and PtN2 (see also Ref. 7) 
or into Pt and Pt 3 N. However, the series of the possible 
phase transitions must be carefully investigated. 



B. Volume per Atom and Lattice Parameters 

The obtained equilibrium volume per atom Vo, i.e. the 
inverse of the number density, for all the considered mod- 
ifications are numerally presented in Table I and graphi- 
cally depicted relative to the Pt(Al) in Fig. 4. On aver- 
age, Pt 3 N phases tend not to change the number density 
of the host parent Pt(Al); PtN phases tend to slightly 
increase it; while the PtN2 increase it significantly. 

It is also evident from Fig. 4 that in crossing the board- 
ers between the Pt 3 N and PtN and between the PtN and 
PtN2 islands, i.e. in increasing the N content, Vb tends to 
decrease while the volume per Pt atom V^ Pt , a measure 
of the average Pt-Pt bond length, tends to increase. The 
latter finding has been found to be true for the nitrides 
of Cu 43 and Ag 61 as well. 



C. Bulk Modulus and its Pressure Derivative 

With only a few exceptions, Fig. 4 and Table I reveal 
that nitridation of Pt apparently tends to reduce its bulk 
modulus. Relative to each other, the twenty B$s show no 
clear trend. The most energetically favorable PtN phase, 
B17, has 42 GPa higher bulk modulus than the proposed 
PtN(B3). 

As we mentioned somewhere else 67 , Bo is far more sen- 
sitive to any change in volume than the change in E co h. 
The case of PtN(B2) is a clear example, in which the 
slight decrease in Vq overcomes the significant increase in 
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E co h leading only to a very small decrease in Bq (Fig. 4 
and Table I). 

Given that all the considered phases have higher -E co h 
than Pt(Al), the foregoing argument fails to explain the 
decrease in Bq in the case of the structures which have 
lower Vq than their parent Pt(Al) and have lower £7 C oh 
than the extreme case PtN(B2). However, if one replaces 
Vq in the argument above with V pt , the contradiction 
can be lifted. Therefore, we believe that the mechanical 
properties in these nitrides may be dominated by the 
effect of the Pt-Pt bond length more than the simple 
number density. 

Although the GGA calculated Bq values in the present 
and previous works (Table I) are far smaller than the 
reported experimental value, our obtained bulk modu- 
lus for PtN 2 (Cl) is 20 GPa higher than that of Pt(Al). 
This is exactly the measured value for Pt after the PtN 
formation took place. The observation was considered 
by Gregoryanz et al. as an indication that some N is 
dissolved in Pt 1 . Recalling that the Bq of the produced 
platinum nitride is ~ 100 GPa than that of Pt(Al) 1 , our 
GGA-obtained B for PtN 2 (Cl) is - 80 GPa less than 
the experimental value 68 . 

It may be worth to notice from Table I that the 
lattice parameter a of PtN2(Cl) is 0.13 A higher than 
that of PtN(B3); yet the B of the former is ~ 70 GPa 
higher than the latter. This difference in Bq can be 
attributed to the fact that in B3, N atoms occupy 
only half of the tetrahcdral interstitial sites of the Pt 
sub-lattice, while in CI, the four remaining tetrahcdral 
interstitial sites are filled with N atoms 8,9,69 . This filling 
significantly reduces the compressibility but slightly 
increases the volume of the unit cell. This fact can 
also be seen readily as a consequence of the difference 
in the average volume per atom in the two cases (Table I). 

The pressure derivative of the bulk modulus, B' , mea- 
sures the sensitivity of Bq to any external pressure. The 
top subfigure in Fig. 4 reveals that the bulk moduli 
of Pt 3 N(Fe 3 N) and PtN 2 (C2, C18 and CoSb 2 ) increase 
upon application of external pressure. Pt3N(D0 2 ) is very 
sensitive and its Bq will increase significantly under an in- 
finitesimal excess of pressure. Pt 3 N(Ll 2 and RhF 3 ) tend 
to be inert; while Pt 3 N(D0 3 , A15 and D0 9 ), PtN 2 (Cl) 
and all PtN phases tend to decrease their bulk modulus 
upon application of external pressure. Although B is a 
measurable quantity 55 , we couldn't find any experimen- 
tal value to test our obtained values against. 



D. Formation Energies 

From Fig. 4 and Table I, it is evident that formation 
energy Ef has the same trend as the cohesive energy 
i? co h- If Ef is taken as a measure of synthesized, then 
the relatively most favorable Pt 3 N phases have the same 
synthesized as the most favorable PtN and PtN 2 . 

A positive value of Ef means, in principle, that, at the 



temperature and pressure at which Ef is calculated, the 
phase is thermodynamically unstable (endothcrmic) and 
have a tendency to decompose into its constituent com- 
ponents. In our case, this observation is corroborated 
by the experimental fact that the synthesis of the plat- 
inum nitrides was achieved only at high temperature and 
temperature 6 ' 60 . 

Using different methods, other researchers 4 ' 6,60 also 
obtained positive (zero-pressure and zero-temperature) 
formation energies for some PtN and/or PtN 2 phases. 
Some of their values are included in Table I with indica- 
tion to the methods of calculations. 

The obtained relative difference in Ef for PtN 2 (Cl) 
and PtN 2 (C2) is in good agreement with Ref. 4. How- 
ever, the differences in our and their obtained Ef values 
can be attributed to three factors: First, the difference 
in the obtained lattice parameter (see Table I). Second, 
the value of our calculated equilibrium free parameter 
u is 0.417 while Ref. 4 obtained 0.415 70 . Third, and 
the most significant source of difference, the experimen- 
tal value of i? co h(N| as ) in Eq. 7 was used by Ref. 4, while 
we calculated it as described in sub-section II D. 

It may be worth mentioning here that a negative the- 
oretical value of Ef = —0.4 eY/atom was obtained for 
PtN 2 (C2) at P = 50 GPa, showing excellent agreement 
with experiment 4 . Moreover, Young et al. 6 claimed that 
PtN 2 dissociates upon mild heating below P = 10 GPa. 



E. Electronic Properties 

The DFT obtained band diagrams ef(k) and spin- 
projected total and partial density of states (DOS) of 
the most stable modifications: Pt 3 N(RhF 3 ), PtN(B3 and 
B17), and PtN 2 (C18) arc displayed in Figs. 6, 7, 8 and 
9, respectively. Spin-projected total density of states 
(TDOS) are shown in sub-figure (b) in each case. Be- 
cause in these four considered cases electrons occupy the 
spin-up and the spin-down bands equally, it was sufficient 
only to display spin-up DOS and spin-up band diagrams. 
Displaying the energy bands along densely sampled high- 
symmetry strings of k-points allows us to extract infor- 
mation about the electronic structure of these phases. 
Moreover, to investigate the details of the orbital char- 
acter of the bands, the Pt(s,p, d) and N(s,p) resolved 
DOS's are plotted at the same energy scale. 

With The Fermi surface crossing the partly occupied 
bands, it is clear from Figs. 6, 7, and 8 that Pt 3 N(RhF 3 ) 
PtN(B3) and PtN(B17) are metals. 

The TDOS of Fig. 9(b) reveal that PtN 2 (C18) is a 
semiconductor with (Fig. 9(a)) its valence band maxi- 
mum (VBM) at (F, —0.091 eV) and its conduction band 
minimum (CBM) at (Y, 0.044 eV), resulting in a narrow 
direct energy band gap E g = 0.135 eV of width. Below 
this fundamental gap there are three bands: the deep one 
at ~ —20.5 eV consists mainly of the N(2s) states. Its 
high DOS and sharp feature correspond to its little and 
slow energy variation in the k— space. The second band 




FIG. 6. (Color online.) DFT calculated electronic structure for Pt3N in the RI1F3 structure: (a) band structure along the 
high-symmetry k-points which are labeled according to Ref. [71]. Their coordinates w.r.t. the reciprocal lattice basis vectors 
are: F(0.5, 0.5, 0.0), Q(0.375, 0.625, 0.0), 5(0.5,0.75,0.25), Z(0.5, 0.5, 0.5), T(0.0, 0.0, 0.0), £(0.0, 0.5, 0.0), F(0.25, 0.5, -.25), 
£(0.0, 0.5, — .5); (b) spin-projected total density of states (TDOS); (c) partial density of states (PDOS) of Pt(s,p, d) orbitals 
in Pt 3 N; and (d) PDOS of N(s,p) orbitals in Pt 3 N. 




FIG. 7. (Color online.) DFT calculated electronic structure for PtN in the B3 structure: (a) band structure along the high- 
symmetry k-points which are labeled according to Ref. [71] a . Their coordinates w.r.t. the reciprocal lattice basis vectors 
are: F(0.0, 0.0, 0.0), X(0.0, 0.5, 0.0), W (0.75, 0.25, 0.5), L(0.5, 0.5, 0.5), A" (0.750, 0.375, 0.375), (7(0.625,0.250,0.625); (b) spin- 
projected total density of states (TDOS); (c) partial density of states (PDOS) of Pt(s,p, d) orbitals in PtN; and (d) PDOS of 
N(s,p) orbitals in PtN. 

a The coordinates of the W point is not as the same as in Ref. [71], but they are equivalent. Also, the coordinates of U and K 
are not given in Ref. [71]. The coordinates of U, K and the equivalent W were created by means of XCrySDen! 




FIG. 8. (Color online.) DFT calculated electronic structure for PtN in the B17 structure: (a) band structure along the high- 
symmetry k-points which are labeled according to Ref. [71]. Their coordinates w.r.t. the reciprocal lattice basis vectors are: 
X(0.0,0.5,0.0), M(0.5,0.5,0.0), T(0.0, 0.0, 0.0), £(0.0,0.0,0.5), ^(0.5,0.5,0.5), R(0.0, 0.5, 0.5); (b) spin-projected total density 
of states (TDOS); (c) partial density of states (PDOS) of Pt(s,p, d) orbitals in PtN; and (d) PDOS of N(s,p) orbitals in PtN. 
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is relatively narrow (~ 2.2 eV of width) with low density 
and steming mainly from a mixture of the N states with 
Pt(cf) states. The superposition Pt(d) and N(p) states 
in the region from —10.314 eV to —0.091 eV below the 
fundamental gap constitutes the third band with highly 
structured, intense and narrow series of peaks. Our ob- 
tained TDOS and PDOS show excellent agreement with 
Ref. 60 where also PtN2(C18) was predicted to be a 
semiconductor, but band diagrams and E g value are not 
given. 

It may be worth mentioning here that PtN(Bl) 10 and 
PtN(B4) 59 were found to be metallic, PtN 2 (Cl) was 
found to be a poor metal 8 , PtN 2 (CoSb 2 ) 60 was found to 
be a semiconductor, and an indirect band gap between 
1.2 eV 60 and 1.5 eV 6 has been obtained for PtN 2 (C2). 



F. Optical Properties 

GW calculations were carried out for the PtN(B3) and 
PtN(B17) metallic phases at their equilibrium. Figs. 10 
and 11 display the obtained real and imaginary parts of 
the frequency-dependent dielectric function £rpa(w) of 
these two phases and the corresponding derived optical 
spectra (Eqs. 12 - 15). In each sub-figure, the optical re- 
gion [~ (3.183 - 1.655) eV = (390 - 750) nm] is shaded. 

The non-vanishing absorption coefficient a (w) in the 
whole range for both phases confirms their metallic char- 
acter. As it should be the case, refraction n(u>) and ex- 
tinction k(lu) coefficients behave as the real e le (ui) and 
the imaginary £i m (w) dielectric functions, respectively. 

As one can see from sub-figure 10(d), close to the edge 
of the optical region at ~ (1.762 eV = 703.768 nm) 
PtN(B3) is 50% reflector and 50% transmitter. From 
- (2.071 eV = 598.579 nm) to the UV region, PtN(B3) is 
only ~ 40% reflecting but ~ 60% transmitting. However, 
more of the transmitted portion in this region will be 
absorbed as the photon energy increases. This fact can 
be readily noticed if one compares sub-figures 10(c) and 
10(d). 

PtN(B17), as can be seen from sub-figure 11(d), 
is a very good reflector in the whole region until 
~ (3.000 eV = 413.281 nm) where it equally reflects and 
transmits the violet light. However, less of the transmit- 
ted portion in the optical region will be absorbed as the 
photon wavelength decreases. This fact can be readily 
observed in sub-figures 11(c). 

According to the best of our knowledge, there is no 
available experimental optical spectra for the platinum 
nitride. However, from their visual appearance, all the 
synthesized platinum nitride samples look very shiny and 
darker than their parent platinum in reflected light and 
totally opaque in transmitted light. These features sug- 
gest that PtN is either a poor metal or a semiconductor 
with a small band gap 1 . 

From Figs. 11 and 8, the above mentioned properties 



are strongly met by PtN(B17), but purely seen (Figs. 10 
and 7) in PtN(B3), as discussed above. Unfortunately, 
we did not carry out optical calculations for PtN 2 (Cl or 
C2). 

G. PtN versus PtN 2 

Using our own obtained results in the present work as 
well as the findings of other researchers, below we make 
a comparison between the PtN modifications (supported 
by the experimentalists) and the PtN 2 phases (supported 
by the theoreticians): 

• Given that GGA calculated lattice parameters are 
usually overestimated 72-74 , the obtained values of 
the a lattice parameter for PtN 2 (Cl and C2) are 
the closest ones to the experimental value (to 
within 3 % and 2 %, respectively), while the PtN 
phases are in poor agreement with experiment, as 
can be seen in Table I. 

• First-principles studies of transition metals nitrides 
show that the Bo's of the elemental metals are 
generally enhanced by nitridation 1 ' . Compared 
to experiment, Table I and Fig. 4 reveal that 
this trend is met by PtN 2 (Cl), while PtN(B3) has 
50 GPa lower than Pt(Al). 

• Like the first synthesized sample and the proposed 
PtN(B3) modification 1 , PtN 2 (Cl 8 and C2 4 ) have 
fee sub-lattice of Pt. 

• PtN 2 ( CI 8 , C2 6 > 60 , C18 60 and CoSb 2 60 ) have 
all been found to be elastically stable, while 
PtN(B3 5 ' 6 ' 8 ' 10 and B17 10 ) were found to be elas- 
tically unstable. 

• Formation and cohesive energies of PtN 2 (C2, C18 
and CoSb 2 ) are lower than that of PtN(B3) [Table 
I and Fig. 4] . 

• In excellent agreement with experiment, the calcu- 
lated formation energy of PtN 2 (C2) at P = 50 GPa 
was calculated to be negative 4 , while calculations 
found PtN(B3) to be thermodynamically unstable 
at all pressures 9 . 

• The experimentally obtained Raman spectrum of 
the reproduced platinum nitride 4 matches closely 
that of pyrite (FeS2), i.e. in the C2 structure, 
but does not match the PtN(B3) spectrum that 
expected from group theory 4 . 

• The theoretically calculated 6 ' 60 Raman spectrum 
for PtN 2 (C2) shows good agreement with the first 
experimentally obtained one 1 . 

• In agreement with the experimental observation 
and the visual appearance of the first produced 
platinum nitrides 1 , PtN 2 (Cl) was found to be a 
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FIG. 9. (Color online.) DFT calculated electronic structure for PtN2 in the C18 structure: (a) band structure along the high- 
symmetry k-points which are labeled according to Ref. [71]. Their coordinates w.r.t. the reciprocal lattice basis vectors are: 
r(0.0, 0.0, 0.0), X(0.0, 0.5, 0.0), 5(-.5, 0.5, 0.0), Y(-.5, 0.0, 0.0), Z(0.0, 0.0, 0.5), U(0.0, 0.5, 0.5), R(-.5, 0.5, 0.5), T(-.5, 0.0, 0.5); 
(b) spin-projected total density of states (TDOS); (c) partial density of states (PDOS) of Pt(s,p, d) orbitals in PtN2; and (d) 
PDOS of N(s,p) orbitals in PtN 2 . 
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FIG. 10. (Color online.) The GW calculated frequency-dependent optical spectra of PtN(B3): (a) the real e Ie (w) and the 
imaginary £i m (w) parts of the dielectric function £rpa(i^); (b) refraction n(cu) and extinction k(lj) coefficients; (c) absorption 
coefficient a(w); and (d) reflectivity R(lu) and transmitivity T(uj). The shaded window highlights the optical region. 



(a) (b) (c) (d) 




FIG. 11. (Color online.) The GW calculated frequency-dependent optical spectra of PtN(B17): (a) the real e ro (oj) and the 
imaginary £i m (w) parts of the dielectric function £rpa(w); (b) refraction n(co) and extinction k(cj) coefficients; (c) absorption 
coefficient a(w); and (d) reflectivity R(lj) and transmitivity T(lj). The shaded window highlights the optical region. 
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poor metal 8 , and we found PtN 2 (C18) to be a semi- 
conductor with a small band gap. 

Hence, in contrast to the proposed PtN modifications, 
PtN 2 phases possess many similar properties as the syn- 
thesized phase 76 . 

IV. CONCLUSIONS 

In summary, we presented a systematic series of first- 
principles calculations of the energy-optimized geome- 
tries, phase stabilities and electronic and optical prop- 
erties of bulk Pt 3 N, PtN and PtN 2 in twenty different 
crystal structures. Comprehensive comparison with ex- 
periment and with previous calculations has been made, 
and excellent agreement has been achieved. We found 
that both the lowest energy and the highest bulk modu- 
lus phases belong to the PtN 2 series and not to the PtN 
family. Moreover, the calculated electronic and optical 
properties of the PtN 2 phases show stronger consistency 
with experiment than the claimed PtN(B3) phase. In 



the present work, we have investigated a wider parame- 
ter sub-space than previous calculations, and to the best 
of our knowledge, the present work is the first to propose 
and to study the physical properties of PtaN, as well 
as the first to theoretically calculate the optical spec- 
tra of this new material. However, optical properties of 
PtN 2 (Cl and C2) have not been investigated, and we 
strongly recommend optical calculations for these phases 
and obtained results should be tested against experiment. 
Moreover, experimentalists should provide the commu- 
nity with more data. 
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